Scattering induced spin orientation and spin currents in gyrotropic structures 
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It is shown that additional contributions both to current-induced spin orientation and to the 
spin Hall effect arise in quantum wells due to gyrotropy of the structures. Microscopically, they 
are related to basic properties of gyrotropic systems, namely, linear in the wave vector terms in 
the matrix element of electron scattering and in the energy spectrum. Calculation shows that in 
high-mobility structures the contribution to the spin Hall current considered here can exceed the 
term originated from the Mott skew scattering. 



PACS numbers: 72.25.Pn, 72.25.Rb, 73.63.Hs 
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I. INTRODUCTION 

Spin-orbit coupling in low-dimensional semiconductor 
structures is attracting a great deal of attention since 
it underlies effects of manipulating spins of charge car- 
riers by electrical means. Spin orientation of free car- 
riers by electric current 0, 0, 0, 13 and the spin Hall 
effect, where a charge current drives a transverse spin 
current [E IS 0> are among the widely discussed phe- 
nomena in the thriving field of semiconductor spintron- 
ics. Microscopically, both effects are caused by spin-orbit 
interaction and can be related either to spin-dependent 
scattering (extrinsic contributions) or solely to spin split- 
ting of the band structure (intrinsic terms) . Comparative 
roles of the extrinsic and intrinsic mechanisms in spin 
transport of charge carriers in bulk and low-dimensional 
semiconductors are, at present, the subject of experimen- 
tal and theoretical discussion (see Refs. 0,0 for review). 

So far, scattering- induced spin effects in transport of 
conduction electrons have been attributed mainly to the 
Mott skew scattering 0, 0, EO, EH The Mott term 
in the matrix element of scattering can be written as 
A er • [k X k'] , where A is a parameter, er is the vector of the 
Pauli matrices, k and k' are the initial and scattered wave 
vectors, respectively. This term is quadratic in the wave 
vector and present in any, even centrosymmetric, struc- 
ture. In semiconductor quantum wells (QWs) the Mott 
term is not dominant contribution to the spin-dependent 
part of electron scattering since gyrotropic symmetry of 
QW structures allows for linear in the wave vector cou- 
pling of spin states. An example of such a fc-linear spin- 
orbit coupling is the Rashba or the Dresselhaus spin-orbit 
splitting of the electron subbands induced by Structure 
and Bulk Inversion Asymmetry (SIA and BIA, respec- 
tively). Similarly, a spin-dependent term linear in the 
wave vector appears in the amplitude of electron scatter- 
ing by static defects or phonons. Taking into considera- 
tion this contribution and neglecting the Mott term, the 
matrix element of electron scattering can be presented as 
(see Ref. 01 an d references therein) 



V k 'k = V + >J V a/3 a a {kp + k'p) , 



(1) 
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where Vq is the matrix element of conventional spin- 



conserving scattering and V a p are parameters determined 
by space distribution and structure details of the scat- 
terers. In the case of elastic scattering from short-range 
static defects, that is assumed below, the parameters V a p 
and Vq are independent of the wave vectors k and k' . 
Linear in the wave vector term seems to be the domi- 
nant spin-dependent contribution to the matrix element 
of electron scattering in QWs. Unlike the Mott contri- 
bution, it can be obtained in first order of the k • p per- 
turbation theory [T^ |. 

In this paper we analyze spin transport of two- 
dimensional (2D) electrons in the presence of fc-linear 
terms in the scattering amplitude. We show that these 
terms together with the spin-orbit spectrum splitting 
give rise to additional contributions to both the current- 
induced spin orientation of free carriers and the spin Hall 
effect. The contribution to the spin orientation of 2D 
electrons is comparable in magnitude to that related to 
current-induced carrier redistribution between the spin- 
split subbands Hi . As regards the spin Hall current, 
the proposed contribution can exceed the term originated 
from the Mott skew scattering in QW structures with 
high mobility. 



II. MICROSCOPIC MODEL 

Microscopically, mechanism of spin Hall current gen- 
eration and spin orientation of 2D electrons by electric 
current due to fc-linear terms in the scattering amplitude 
is as follows. Application of an electric field E in the 
QW plane results in a directed flow of the carriers. Due to 
spin-dependent asymmetry of scattering, electrons driven 
by the electric field are scattered in preferred directions 
depending on their spin states. It leads to generation of 
a spin-dependent electron distribution, where particles 
with a certain wave vector fc carry a certain spin ori- 
entation. The explicit form of the distribution depends 
on the origin of the spin-dependent scattering. Since fc- 
linear terms in the scattering amplitude are caused by 
inversion asymmetry of QWs, one can distinguish, simi- 
larly to the spectrum splitting, the SIA and BIA contri- 
butions to the scattering amplitude. The corresponding 
distributions of the spin density in fc space are shown in 
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Fig. 1 for electrons confined in (OOl)-grown QW and sub- 
jected to the in-plane electric field along [f 10] and [110] 
axes. Detailed calculation of the presented distributions 
is given in the next section. Here we only note that the 
angular dependencies of the spin distributions are de- 
scribed by quadratic harmonics. Average spin polariza- 
tion remains zero and, in contrast to the Mott scattering, 
electrons with the opposite wave vectors k and —k carry 
the same spin. Therefore, the spin-dependent scattering 
of the form does not lead itself to a current-induced 
spin orientation nor to the spin Hall effect. 

A net spin orientation of the electron gas and a spin 
current appear as a result of the subsequent spin dynam- 
ics of carriers. The spin dynamics of the conduction elec- 
trons is known to be governed by spin-orbit coupling that 
may be considered as an effective magnetic field acting 
on electron spins. In this effective magnetic field spins 
of the carriers, initially directed according to the spin- 
dependent scattering processes, precess. Such a spin dy- 
namics is illustrated in Fig. 2. To be specific we consider 
that the electric field E is directed along the [110] axis 
and both the effective magnetic field and the fe-linear 
terms in the scattering amplitude are induced by struc- 
ture inversion asymmetry. The effective magnetic field 
caused by spin-orbit coupling in QWs is known to be an 
odd function of the wave vector k. In the the case of SIA 
it is directed perpendicular to k as shown in Fig. 2 by 
dashed arrows. Therefore, electrons moving along or op- 
posite to the x|| [110] axis carry spins oriented parallel (or 
antiparallel) to the effective magnetic field, while the par- 
ticles moving along or opposite to the j/ 1 1 [110] axis carry 
spins oriented perpendicular to the effective field. As a 




FIG. 1: Distribution of electron spin caused by fe-linear terms 
in the scattering amplitude for various directions of the elec- 
tric field. Figs. (a),(b) and (c),(d) correspond to the SIA and 
BIA fe-linear terms in the scattering amplitude in (OOl)-grown 
QWs, respectively. 
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FIG. 2: Microscopic mechanism of the current-induced spin 
orientation and the spin Hall effect. Spin-dependent asym- 
metry of electron scattering followed by spin precession in 
the effective magnetic field leads to (i) generation of a spin 
current and (ii) net spin orientation of carriers. 



result of the spin precession, the spin component S z > 
(S z < 0) appears for the carriers with the positive (neg- 
ative) wave vector k y . This state corresponds to the spin 
Hall effect where an electric field drives a transverse spin 
current, i.e. oppositely directed flows of carriers with the 
opposite spins. Moreover, the spin precession in the ef- 
fective magnetic field leads to appearance of an average 
spin orientation of carriers in the QW plane, opposite 
to the y axis. Indeed, electrons with the spins directed 
opposite to the y axis are affected by collinear magnetic 
field and retain the spin orientation, while the carriers 
with the spins directed along the y axis partially lose the 
polarization due to the spin precession. The rate of the 
spin generation is determined by the average angle of spin 
rotation in the effective magnetic field, similarly to the 
Hanle effect. 



III. THEORY 

Theory of the scattering-induced spin orientation and 
the spin Hall effect is developed here by using the spin- 
density-matrix technique. Dynamics of the density ma- 
trix pk of electrons subjected to an in-plane electric field 
E is given by the kinetic equation 

dpk , „dp k i 

-dT + eE dhk + h [Hs °> Pk] = Stpk - (2) 

Here e is the electron charge, H so is the Hamiltonian of 
spin-orbit coupling that describes spin precession in an 
effective magnetic field 

H so = ^2 7<*0 Vakp = - Ofe • <T , (3) 

a/3 
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"f a /3 are the structure parameters, fife is the Larmor fre- 
quency corresponding to the effective field, and Stpk is 
the collision integral. In the present letter we restrict 
ourselves to the case of low temperatures when transport 
properties are determined by the Fermi surface and en- 
ergy mixing of carriers is suppressed. For elastic scatter- 
ing from static defects uniformly distributed in the QW 
plane the collision integral has the form 

Stpk = ■rNd'^2(2Vkk' Pk'V k ' k - Vkk'Vk'kPk (4) 

k' 

-RkVkk'Vk'k) S(e k - Sk') , 

where Nd is the sheet density of defects. Note, that spin- 
orbit splitting of the energy spectrum is neglected in the 
collision integral J3J, since the splitting is much smaller 
than the electron kinetic energy. Corrections to the en- 
ergy due to the spin-orbit splitting as well as similar cor- 
rections to the electron velocity, which can be crucial for 
intrinsic spin effects, are unimportant here. 

The density matrix can be presented as follows 



Pk = (/o + Sfk) I + Sk ■ a- 



(5) 



where fo is the function of the equilibrium carrier dis- 
tribution, / is the 2x2 unit matrix, Sfk and Sk are 
the electric field-induced corrections to the diagonal and 
spin components of the density matrix. We assume that 
the electric field oscillates at the frequency uj, E cx 
exp(— iojt). Then, in linear in the electric field regime, 
the terms Sfk and Sk have the same time dependence. 
To first order in the parameters V a p the charge transport 
is independent of the spin part of scattering amplitude 
and the diagonal correction to the equilibrium density 
matrix has the form 



Sfk 



er E ■ v dfo 
1 — iojt de 



(6) 



where v = Tik/m* is the velocity, m* is the electron 
effective mass, r = h 3 /(V 2 Ndm*) is the momentum 
isotropization time, and e is the electron kinetic energy. 
Note, that Sfk is the correction that describes the con- 
ventional (Drude) ac conductivity of the electron gas. 

Equation for the spin component of the density matrix 
Sk can be derived from Eq. J2J). To first order in spin- 
dependent part of the matrix element of scattering V s it 
takes the form 

Sk — Sk 

- iujSk + [S k x fl k ] = h g k , (7) 



collisions with structure defects, while the right-hand side 
of the equation stands for the scattering-induced spin re- 
distribution among electron states with different wave 
vectors. Particularly, the first term on the right-hand 
side of Eq. J7J) describes isotropization of the spin den- 
sity and the second term corresponds to the spin den- 
sity generation due to carrier drift in the electric field in 
the presence of spin-dependent scattering. We note, that 
spin-flip processes are neglected in Eq. iJSJ since they ap- 
pear only in second order in V s . For linear in the wave 
vector terms in the matrix element of scattering given by 
Eq. JQ) , components of the spin generation rate have the 
form 



9k,c 



2eh/m* dfo 



1 



'mot de 



V^_ 
V 



k 2 



E„ 



(9) 



Dependence of the spin generation gk on the wave vec- 
tor is determined by both the direction of the applied 
electric field E and the explicit form of the coefficients 
V a f). The latter is governed by the QW symmetry and 
can be varied. Symmetry analysis shows that in (001)- 
grown QWs there are two nonzero components of the 
pseudotensor V a p, namely V xy and V yx , which can be 
expressed via the SI A and BlA contributions as follows: 
V xy = Vbia+Vsia, V yx = Vbia-Vsia- Figure 1 demon- 
strates scattering-induced distributions of the spin polar- 
ization gk in k space for the electric field directed along 
the x and y axes. Figures (a),(b) and (c),(d) correspond 
to the cases where the fc-linear terms in the matrix ele- 
ment of scattering are caused by SI A (V xy = — V yx ) and 
BIA (V xy = V yx ), respectively. 

As mention above, the net spin orientation of carriers 
and the spin Hall current appear as a result of the spin- 
dependent scattering and the subsequent precession of 
electron spins in the effective magnetic field. Following 
Eq. (J7J and taking into account that the frequency f2fc is 
an odd function of the wave vector, while gk is an even 
function of k and gk — 0, one can derive the equation for 
the spin component Sk 



ft k x [Ofc x (Sk + rgk)} 
[(1 - iwr) 2 + (Q fe r) 2 ] 



iuSk = . 



(10) 



where the angle brackets mean averaging over directions 
of the wave vector. Solutions of Eqs. (7J and (|10fl allow 
one to find the spin components of the density matrix 
and calculate the spin polarization of the electron gas as 
well as the spin Hall current. 



where Sk is Sk averaged over directions of the wave vec- 
tor k, gk is the spin generation rate into the state with 
the wave vector k, 

9k = ^-N d Tr[<rV V s ](6f k , - Sf k ) S(e k - e y ) . (8) 
n k 

The left-hand side of Eq. J7J describes spin dynamics of 
the carriers during their free motion between consecutive 



IV. RESULTS AND DISCUSSION 

Calculation shows that in the case Sl^r <C 1 the aver- 
age electron spin, defined as J2k Sk/N e , with N e being 
the carrier concentration, has the form 



s x (w) 



V l yx 



V 



K V (U) 



1 - \wT x (l 



(11) 
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Vxy 'fyx ^7 
V Jxy 



y% 
V 



K x (oj) 



1 — iujT y (l — ujt) s 



Here V xy , V yx , j xy and "f yx are nonzero components of 
the pseudotensors V a p and in (OOl)-grown structures, 
K(lo) is the average electron wave vector corresponding 
to the drift velocity of carriers in the electric field 



er/h 



1 



(12) 



and T a (a = x,y) are the D'yakonov-Perel' relaxation 
times of the spin components 



l/T a = r[(n|)-(fi? }] 



(13) 



Magnitude of the spin orientation achieved by dc electric 
current s(Q) depends on the ratio J x y/jyx rather than on 
absolute value of the spin-orbit splitting Ml^ . This can 
be understood considering the average spin as a balance 
between processes of spin generation and spin relaxation. 
s(0) is given by Ts, where the spin relaxation time T ex 
and the spin generation rate s oc f2|. 
The spin current is characterized by a pseudotensor 
J with the components Jp describing the flow in the (3 
direction of spins oriented along the a axis. In terms of 
the kinetic theory such a component of the spin current 
is contributed by a non-equilibrium correction oc a a kp to 
the electron spin density matrix and given by 



V{3 



(14) 



Calculation shows that the components of the scattering- 
induced spin Hall current in (OOl)-grown QWs have the 
form 



4» 
4» 



' 2% 
1 2% 



Vyx 



V 

Vyx 

Vn 



Vxv 



ly* 



V 

Vxy 
Vn 



ojT x j y (oj)/e 



7yx 



— iu)T x (l — VjJt) 
ojT y j x (oj)/e 
1 — iujTyil — iojt) 2 



where hp is the Fermi wave vector and j{oj) = 
eN e hK (oj) / m* is the charge current. 

Equations (|llfl and (|15|l describe the current-induced 
spin orientation and the spin Hall current for the case 
fifeT -C 1. In high-mobility 2D structures this inequal- 
ity may be violated. At arbitrary fifcr solutions of 
Eqs. (0) and (|10|) have complicated form but become 
simpler if the absolute value of the Larmor frequency 
fife is independent of the wave vector direction. This 
is fulfilled when the spin-orbit splitting is caused only 
by BIA or by SIA, "f xy /"f yx = ±1, respectively. In 
this particular case the spin orientation and the spin 
Hall current are given by Eqs. (|ll|l and l)15|l. where 
the denominators 1 — iuT a (l — \ujt) 2 are replaced by 
1 - iwT[(l - iojt) 2 + 2t/T], with the time T (T x = T y ) 
being determined by Eq. (|13J) . 

Frequency dependencies of the current-induced 
spin polarization |s y (w)| and the spin Hall current 
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FIG. 3: Frequency dependencies of (a) the spin orientation 
|sj,(a;)| and (b) the spin current \eJy(uj)/j m (0)\ for different 



Qk F r and V xv /Vo = 10' 
10 5 cm" 1 . 



k F = 10 b 



K(0) = 



\eJy(u))/j x (0)\ are plotted in Fig. 3 for different pa- 
rameters f2fc F T. It is assumed that the electric field is 
directed along the x axis and the fe-linear terms both in 
the scattering amplitude and in the spectrum splitting 
are related to structure inversion asymmetry. Magnitude 
of the spin orientation reaches maximum in the dc limit 
and decreases with increasing the field frequency u>, see 
Fig. 3(a). In systems where fi^r <C 1, the spin relaxation 
time is much longer than the momentum relaxation 
time, T 3> r, and the frequency dependence of the 
spin orientation s a (uj) is given by s Q (0)/(l — iwT Q ). 
Such a behavior is natural because the spin orientation 
represents a carrier redistribution between the spin 
states and the spin dynamics is governed by the spin 
relaxation time. 

In contrast to the current-induced spin orientation, the 
spin Hall current reveals nonmonotonic frequency depen- 
dence. It vanishes at zero frequency, increases linearly 
with the frequency at low lj, reaches maximum and de- 
creases when oj exceeds 1/r. The cancellation of the spin 
Hall effect in the dc limit as well as the linear growth 
at small oj is a direct consequence of the particular re- 
lation between spin dynamics and spin fluxes in systems 
with spin-orbit splitting of the band structure linear in 
the wave vector |l3j . The relation between components 
of the spin and the spin current can be easily obtained 
from Eq. J7J. Summing Eq. (J7J over k and taking into 
account that f2fc is a linear function of the wave vector, 
in particular, Qk^ x oc k y and £lk, y oc k x in (OOl)-grown 
structures, one obtains the relation J^(u>) oc ojs a (oj). At 
small frequencies oj the spin orientation is a finite value, 
implying the linear growth of the spin Hall current with 
the frequency. In the range l/T <C oj <C 1/t the spin 
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orientation drops asw" 1 and the frequency dependence 
of the spin Hall current has a plateau. After that, in 
the high-frequency limit, lu 3> 1/t, the spin orientation 
and the spin Hall current decrease as w -4 and uj~ 3 , re- 
spectively. We note that in the high-frequency range the 
absorption of ac electric field by free carriers can be sig- 
nificant and result in pure spin photocurrents [T3.ll5| and 
spin orientation of the carriers |12j . 

Finally, we present estimations for the considered ef- 
fects. Following Eq. {TJ| the current-induced spin ori- 
entation at zero frequency can be estimated as for 
V xy /Vo ~ 1CP 8 cm, appropriate to GaAs-based struc- 
tures, and the drift wave vector K(0) ~ 10 5 cm -1 . This 
value corresponds to the spin orientation caused by other 
mechanism, namely, current-induced electron redistribu- 
tion between the spin-split subbands. For this particu- 



lar mechanism the spin orientation can be estimated as 
jxyKiO)/ Ep that also gives 1CT 3 for j xy ~ 1(T 7 meV •cm 
and the Fermi energy Ep ~ 10 meV [ifj. Further in- 
crease of the Fermi energy (or temperature for the Boltz- 
mann statistics) leads to predominance of the scattering- 
related mechanism. As regards the spin Hall effect, 
the ratio of the spin current (|15(l to the contribution 
caused by the Mott skew scattering can be estimated 
as V xy /{\kp) flk F T. In high mobility structures, where 
£lk F T ~ 1) the spin Hall current caused by fc-linear terms 
in the scattering amplitude exceeds by an order of mag- 
nitude the contribution related to the skew scattering. 
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